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ABSTRACT 

We address the issue of ambipolar diffusion in Herbig-Haro jets. The current 
consensus holds that these jets are launched and collimated via MHD forces. 
Observations have, however, shown that the jets can be mildly to weakly ionized. 
Beginning with a simple model for cylindrical equilibrium between neutral, plasma and 
magnetic pressures we calculate the characteristic time-scale for ambipolar diffusion. 
Our results show that a significant fraction of HH jets will have ambipolar diffusion 
time-scales equivalent to, or less than the dynamical time-scales. This implies that 
MHD equilibria established at the base of a HH jet may not be maintained as the jet 
propagates far from its source. For typical jet parameters one finds that the length 
scale where ambipolar diffusion should become significant corresponds to the typical 
size of large (parsec) scale jets. We discuss the significance of these results for the issue 
of magnetic fields in parsec-scale jets. 



Subject headings: ISM: jets and outflows — magnetic fields — magnetohydrodynamics: 
MHD 



1. INTRODUCTION 

Narrow hypersonic jets are a ubiquitous phenomena associated with star formation. In 
spite of the large database of multi-wavelength observations and numerous theoretical studies, a 
number of fundamental questions remain unanswered about these jets. Paramount among the 
outstanding issues is the role of magnetic fields in the launching, collimation and propagation of 
Young Stellar Object (YSO) jets. The current consensus holds that the jets are formed on small 



scales ( L < 10 AU) via magneto-centrifugal processes associated with accretion disks ( Burrows ei 



al 1996| , |Shu et al. 1994] , Puyed fc Pudritz 1997] ). At larger scales, where jet propagation rather 



than collimation is the issue, there is an implicit assumption that the magnetic fields involved in 
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the launching process will remain embedded in the jets as they traverse the inter-cloud medium. 
Observations of narrowly collimated jets, or chains of bow-shocks, extending out to parsec scale 
distances ( [Reipurth et al. 1997 ) has strengthened the viewpoint that the jet beams must "carry 



their own collimators". This is based on the possibility that without dynamically significant 
magnetic fields confining the beam, the jets may be disrupted by hydrodynamical instabilities 



QStone et al. 1997j , Pardee et al. 1997|) . 



Do the magnetic fields remain embedded in the jets? That is the question we address in 
this paper. If we accept that jets are created via MHD forces then the only way to lose the 
imposed fields is via reconnection or through diffusive processes. The helical topology expected 
for magneto-centrifugally launched jets would not be likely to lead to large scale reconnection 
throughout the beam (though reconnection may have important effects on the radiative properties 
of the jets, Gardiner et al. 1999| ). Thus, if magnetic fields can be cleared out of a jet this is more 



likely to occur via diffusive processes. Large format CCD mosaics have recently demonstrated 
that YSO jets can extend over multi-parsec length scales (Lj « 3 pc). The dynamical age 
(tdyn = Lj/Vj) for these jets falls in the range td yn = 10 4 — 10 5 y ( [Reipurth et al. 1997| , Eisoffcl 



|fc Mundt 1997[ ). Thus, even if diffusive processes are slow, they will have a relatively long time to 
affect the propagation of the beams. In plasmas with less than full ionization, ambipolar diffusion 
is one of the most effective means of driving the rearrangement of flux relative to the plasma. 
Observations of optical Herbig-Haro (HH) jets have shown them to be moderately ionized, with 



ionization fractions ranging from x = .5 — .01 ( Hartigan et al. 1996 , Baccioti fc Eis5ffel 1999 ) 



where the ionization fraction is defined as x = rii/n n and n n and rii are the number densities of 
neutral and ionized particles respectively. The ionization fractions are seen to fall as the plasma 
propagates away from internal shocks and may, therefore, fall to even lower values in the optically 
invisible regions of the beam. Given the potentially low ionization fraction and the extreme ratio 
of jet radius to length, Rj/Lj < 10~ 3 , it may be possible for ambipolar diffusion to significantly 
alter the MHD equilibria established when the jet was launched. 

In this paper we wish to address the issue of ambipolar diffusion in YSO jets. Our present goal 
is simply to establish the characteristic time and length scales for ambipolar diffusion to operate 
with respect to the fundamental jet parameters. We leave its consequences to later papers. In 
section 2 we establish the conditions for an MHD equilibrium in a simplified model of a jet in terms 
of plasma, neutral and magnetic pressures. In section 3 we consider the equation for ion-neutral 
drift in this context and in section 4 we present an equation for the ambipolar diffusion time-scale 
in YSO jets. In section 5 we discuss our results in light of observations and other models. 



2. Initial Configuration 

We begin by considering the simplest model of a magnetized jet. We imagine a cylindrical 
column of material of radius Rj. The column is composed of ions, electrons and neutral species 
and is initially held together internally by a radial balance of pressure (Pi, P e , P n ) and magnetic 
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(J x B) forces. We further assume that the column is in force balance with a (possibly magnetized) 
ambient medium. Thus, initially, the jet does not expand laterally {i.e., in the radial direction). To 
simplify our calculation we assume that the temperature of all species is the same (Tj =T e = T n ). 
We further assume that the magnetic field in the jet is purely poloidal, B = B z {r)e z . Figure 1 
shows a schematic of our model for the magnetized jet. 



2.1. Global Force Balance 

Our calculation is similar to derivations of the ambipolar diffusion time-scale in collapsing 
molecular clouds (a description of those models can be found in Spitzer 1978 and [Mouschovias 



1991]). Molecular clouds are cold gravitationally bound systems of extremely low ionization 
(x ~ 10~ 6 ). This implies that terms involving gas pressure and ion density can be dropped from 
consideration when calculating the effects of ambipolar diffusion. Gravity does not, however, play 
a role in maintaining the cross-sectional properties of a YSO jet and pressure forces can not be 
ignored. In addition, given the range of ionization fractions in these systems the ion inertia should 
not, a priori, be dropped from calculations. In what follows we derive a simple expression for the 
ambipolar diffusion time-scale (based on the ion-neutral drift velocity) including the effects of gas 
and plasma pressure as well as the ion density. 

We begin with the three fluid force balance equations for electrons (e), ions (i) and neutrals 
(n) which can be written as 

(111 —* f — * 1) — * \ —> — * 

Pe-j^ = -VP e -n e e(E + ^xB)+F ei + F eni (1) 



dt \ c 

Pi lf = -V P i+ n i Ze (E + ^ X ^)+Fie + Fin, (2) 

Pn-JT = -VP„ + F„ e + F™, (3) 
at 

where the term d/dt represents the convective derivative and F\2 represents the drag force 
imparted on species 1 by colliding with species 2. Note that F\2 = — ^21 ■ 

In what follows we will ignore the electron inertia term. This is valid when the time-scales 
of interest are long compared with the response time of the electrons. Formally this means that 
the dynamical time t& yn is longer than either the electron cyclotron period 27r/f2 ec or the electron 
plasma period 2-K/u) ep ( Elliott 1993| ) as is the case in YSO jets with fields greater than the \iG level. 



We will also neglect the electron collisional coupling to the the neutrals, F en , given the difference 
in masses, M e « Mi,M n (see footnote #4 Mouschovias 1996| , Ciolek private communication 



1999). 

With these caveats the addition of equations 1 and 2 yields the following expression: 

dvi - - friiZevi- n e ev e B \ - 
Pi — = -V(P< + P e ) + {riiZe - n e e)E + I x B\ + F in . (4) 
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Assuming charge neutrality (riiZ = n e ) yields the following form for the current density, 

J = riiZevi - n e ev e = n e e(vi - v e ), (5) 

and the term with the electric field drops out of equation Making use of Ampere's law, the 
Lorentz force can be decomposed into two terms, 

-(JxB) = -Vi — + -r~^n (6) 

where h is a unit vector directed toward the local center of curvature of the field line, R c is the 
local radius of curvature, and Vj_ is the gradient normal to field lines. 

Our choice of a B = B z (r)e z implies that the second term in equation || is zero. We are left, 
therefore, with only the first term which can be identified as the gradient of magnetic pressure, 
P B = B 2 /8tt. 

Finally, if we add the momentum equations for all three species, again ignoring the electron 
inertia we arrive at the following equation, 

Pi^-+Pn^ = -V(P t + P e + P n + P B ). (7) 

In what follows we need only consider the radial component of the force equation. Thus if, 

^{Pi + P e + P n + P B ) = 0, (8) 

and Vi r = v n ^ r = v e ^ r = (where the v r 's refer to radial velocities) at t = then the entire 
configuration begins in equilibrium. It is, however, what we shall call a quasi- equilibrium since 
each species will be accelerated separately by the terms on the RHS of equations 1-3. The plasma 
and the neutrals will attempt to move past each other as each responds to its own forces. It is 
only the collisional drag Fij between species which holds the initial configuration together. This is 
exactly the situation which occurs in a molecular cloud where the cloud can be supported for some 
period of time by its magnetic field. As in the molecular cloud case, the quasi-equilibrium in the 
jet can not persist indefinitely and the initial configuration will change on a time-scale t a d = Rj/vd 
where Vd is the ion- neutral drift velocity Vd = v n — Vi. We seek to calculate this time-scale. 

Note that equation |8| implies that the pressure distributions of the ions and neutrals is not 
independent. 

§- r Pn = -^ r (P l + Pe + PB) (9) 



We will use this relation below in our calculation of t a d- 
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3. The Ion-Neutral Drift 

If we divide equation |] by pi and subtract from it equation 3 (itself divided by p n ) we find, 

§ - ^ = -iv ( p, + p e + p b ) + -ivp n + P m (E^i) . (10) 

at dt pi p n \ p n pi ) 

If we now assume that the collisions will quickly bring the ion-neutral drift to a steady velocity, 
i.e., the relative accelerations approach zero after a few collisions, we have the following radial 
balance equation between collisions and hydro-magnetic forces, 

F m { !±±»\ = ±|-(J* + P e + Ps) - ^P n . (11) 
V PnPi J Pi or p n dr 

Substituting equation || into the above expression we find the densities cancel, leading to 

F in = ^- r {Pi+Pe + P B ). (12) 



We now make the approximation that the scale of the gradients in the jet are such that 

dP k _ P k 
dr ~ Rj 



(13) 



where Pk is is the characteristic scale for that pressure of the k-th component of the configuration. 
We also use the definition of the plasma (5 parameter, (3 = Pg/Ps, where P g = Pi + P e . Thus 

F in = ^09+1) (14) 

To reduce this equation further we must consider the form of the collisional force. Fi n can be 
written as, 

Fm = PWn M^W n {Vn ~ Vl) ' (16) 

where (crw)in ~3x 10~ 9 cm 3 s _1 is the average collision rate ( |5pitzer 197§| ). We can now solve 
for the ion-neutral drift velocity. 

Mj + M n 1 ( B^_ \ 
Vd = = t r— — (/3 + 1). 17) 

87T {(JW) in \Pip n Rj J 

Thus the time-scale for changes in the initial configuration can be written as 

87r(aw) m M n M (n n RA 2 t 
tad ~ M t + M n X {~^~) ' (18) 

where we have used x = rii/n n . The above expression gives the ambipolar diffusion time-scale in 
terms of the fundamental parameters in the jet (x,n n ,Rj and B). 
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4. Ambipolar Diffusion in YSO Jets 



Of the four variables in equation 18 only first three (x 



n n , Rj) are well characterized 
for HH jets. From observations typical values are: •01 < x < .1; 10 3 < n n j cm < 10 ; 
1 x 10 15 < Rj/ cm < 5 X 10 15 (|Baccioti fc Eisoffel 19991) . Magnetic fields however are not so 



easily categorized. To date there have been only a handful of measurements of fields in YSO jets 
( Ray et al. 1997| ). Thus, it would be better to cast equation |l^ in terms of a parameter such as 
the temperature which has been measured in many YSO jets. Typical values in optical jets are 
5 X 10 3 if < T < 3 X 10 4 K, QBaccioti fc Eisoffel 1999|) . To express equation [l| in terms of T we 



use our assumption that the electron and ion temperature is the same and, furthermore, assume 
that the plasma is mainly composed of Hydrogen atoms. Thus 



Mi = M n 
m = n e 
Pp = n P kT = riikT = Pi 



(19) 
(20) 
(21) 



The magnetic field can now be expressed as 



B 2 = 8vr 



Pe + Pi 

P 



and the ambipolar diffusion time-scale becomes 

tad 



WnriikT—, 



M n {aw) in ( n n R 2 \ ( (3 



4A: 



T; 



(3 + 1 



(22) 



(23) 



Upon normalization the ambipolar diffusion time-scale can be written as 



tad 

Q(P) 



28, 904 



P 
(3+1 



n, 



10 3 cm~ 3 J VlO 15 



( 



'cm 



Q((3)y 



(24) 
(25) 



Note that Q{(3) is positive with an asymptotic value of 1 as (3 — > oo. 



This shows that the ion-neutral drift can be driven by gas pressure forces alone if the 
individual species' pressure gradients have different signs, an artificial situation not likely to be 
encountered in real jets. Detailed calculations of equilibrium magnetically confined jets arising 
from accretion disks indicate that (3 can become quite low in some regions of the jet. As an 
example, we provide in Figure 2 a plot of (3 versus radius for an equilibrium MHD jet. This model 



was calculated via the Given Inner Geometry method of Lery et al. ( Lery et al. 1998 , Lery ei 



al. 1999) for determining the asymptotic structure of magneto-centrifugally launched flows. Figure 
2 shows that (3 can drop to values below = 1 in the inner regions of the jet (r < Rj). 
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In Figure 3 and 4 we show surface plots of t a d vs (5 for = 1 and (3 = .1. From equations 25 
and 26 and these figures it is clear that depending on the conditions in the jets, t a d can become 
as large as 10 6 y and as small as 10 3 y. For jet parameters in the middle of the expected range of 
variation we find t a d of order 10 5 to 10 4 y. What is noteworthy about these results is a significant 
region of parameter space exists where the dynamical time-scale for YSO jets is of order of, or 
greater than, the ambipolar diffusion time: td yn > tad- Thus ambipolar diffusion is likely to play 
a role in the dynamics of large scale YSO jets and outflows. We discuss the consequences of this 
conclusion in the next section. 



5. Discussion and Conclusions 

Our conclusion that ambipolar diffusion time-scales can be comparable to jet dynamical 
lifetimes raises a number of intriguing issues. The first is the most obvious. Ambipolar diffusion 
will rearrange the mass to flux ratio in the jet and alter any initial equilibrium between the 
plasma, neutrals and magnetic field. Our analysis is too simplified to yield conclusions about how 
the jet will evolve in the presence of ambipolar diffusion. In our analysis we did not consider the 
effect of toroidal fields. The hoop stresses associated with toroidal fields will the pull ions towards 
the jet axis. Thus depending on the orientation of the magnetic pressure gradients, magnetic 
forces (pressure and tension) can either compete or apply forces in the same direction. Including a 
toroidal component does not, in general, change the order of magnitude of the ambipolar diffusion 
time-scale, however, it can change the direction in which the ions are pushed. The plasma and 
magnetic field may bleed out of the jet leaving only the neutrals or, conversely, strong hoop 
stresses could draw the plasma and field in towards the jet axis. In both cases, however, the 
neutrals will be free to contract or expand depending on their pressure distribution relative to 
the ambient medium. Thus one expects a potentially significant rearrangement of the jet's cross 
sectional properties when ambipolar and jet time-scales are comparable. The consequences of 
ambipolar diffusion on long term jet propagation should, therefore, be investigated in detail. 

The possibility that many jets will have td yn ~ t a d is suggestive for the general issue of YSO 
outflows. There remains considerable debate over the connection between HH jets and molecular 
outflows ( Masson fc Chernin 1993] , Cabrit 1997| ). The discovery of parsec scale "superjets" 



( Reipurth et al. 1997 , Eisoffel &, Mundt 1997] ) has strengthened the case for jets as the source 



of molecular outflows by extending jet lifetimes. Still, most outflows have larger dynamical 
timescales than even the superjets by a factor of at least a few. One may wonder then why the 
jets have shorter lifetimes then the molecular outflows. Our results suggest that in some systems 
at least it is ambipolar diffusion which sets the lifetime of the visible collimated jet. Converting 
to consideration of length scales, if one takes the minimum ambipolar diffusion time from our 
calculation, t a d ~ 10 3 y and the minimum characteristic jet velocity, Vj ~ 10 2 km/s one finds the 
minimum distance where ambipolar diffusion becomes effective: D a d ~ 0.3 pc Our results indicate 
that ambipolar diffusion should not be effective before a jet reaches this distance. Moreover for 



-8- 



typical values t a d ~ 10 4 y and Vj ~ 3 x 10 2 km/s one finds D a d « 3 pc. It is indeed noteworthy 
that this value corresponds well with the distance where the longest jets start to fade away. When 
tdyn becomes comparable to t a d the beam may no longer be confined by magnetic forces. 

We wish to thank Glenn Ciolek, Larry Heifer and Al Simon for their help with this project. 
We also acknowledge the helpful comments of Tom Ray, Francesca Bacciotti and Alex Raga. This 
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Fig. 1. — Cartoon of Jet Model. Shown is a section of a magnetized jet composed of a plasma, 
neutral gas and magnetic field oriented along the axis of the jet. The pressure and field distributions 
are a function of cylindrical radius only. 
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Beta = 1 . 




2 10 



Fig. 3. — The Ambipolar diffusion time-scale as a function of temperature (in units of 10,000 K) 
and density (in units of 1000 cm -3 ) when [3 = 1. The jet radius is 2 x 10 15 cm. 
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Beta = .1 
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Fig. 4. — The Ambipolar diffusion time-scale as a function of temperature (in units of 10,000 K) 
and density (in units of 1000 cm -3 ) when [3 = .1. The jet radius is 2 x 10 15 cm. 
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